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abstract 

We study stability of noncommutative spaces in matrix models and discuss the contin- 
uum limit which leads to noncommutative Yang-Mills theories (NCYM). It turns out that 
most of noncommutative spaces in bosonic models are unstable. This indicates perturba- 
tive instability of fuzzy MP pointed out by Van Raamsdonk and Armoni et al. persists 
to nonperturbative level in these cases. In this sense, these bosonic NCYM are not well- 
defined, or at least their matrix model formulations studied in this paper do not work. We 
also show that noncommutative backgrounds are stable in a supersymmetric matrix model 
deformed by a cubic Myers term, though the deformation itself breaks supersymmetry. 
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1 Introduction 



Yang-Mills theory on a noncommutative space (noncommutative Yang-Mills theory, 
or simply NCYM) has attracted much interest in theoretical physics. It appears as an 
effective theory of string theory or its matrix models around certain flux backgrounds 
[I] [21 [3] [H m [6]. NCYM contains some interesting physical properties like spacetime 
uncertainty and peculiar solitonic solutions [7]. We also notice that it naturally contains 
gravity (for recent progress, see e.g. [8], [9]). To understand the nonperturbative aspects 
of NCYM better, we need the nonperturbative formulation of it. Matrix models are 
expected to be the most promising approach. Using a matrix model, NCYM is realized 
as an effective theory of a matrix model around a certain background. However, such 
backgrounds are unstable for some cases and whether the theories are well defined or not 
is a nontrivial question. It is well defined only when the backgrounds are stable. In this 
note, we will discuss stability of noncommutative spaces and argue what kinds of NCYM 
can be realized using matrix models. 

Realization of NCYM in matrix models is of interest also from emergent geometry 
point of view. The origin of this concept goes back to early 1980s. The first example, as 
far as we know, is large- reduction [TOl [TT| [T2] which claims that large- gauge theories 
are equivalent to their one point reduced models. In these models, spacetime is embedded 
in gauge fields [HI [121 [13] • We can also find it in the context of quantum theory of 
gravity. From this point of view, spacetime should emerge as a result of some dynamical 
mechanism. As nonperturbative formulations of string theory, various matrix models are 
proposed [2], [3] and, especially in IKKT matrix model [3] , various interpretations are given 
to reahze emergent geometry [TU [151 M- This concept is also discussed in the context of 
AdS/CFT [16] [IT]. 

NCYM is another example of emergent geometry. Let us briefly explain how it shows 
up and what kind of double scaling limit is necessary. We only consider NCYM on a fiat 
noncommutative space and mainly take the continuum limit in which the noncommuta- 
tivity parameter 6 is fixed. We set the gauge group to be U{1) unless otherwise mentioned 
but generalization to U (n) is straightforward. 

For concreteness, let us consider zero-dimensional SU{N) matrix models with a twisted 
boundary condition [12] or a Myers term added [18]. For these models, it is known that 
compact noncommutative spaces like fuzzy spheres are classical solutions. Once we fix 
6, volume of the space and the UV cutoff are related to the matrix size A^. Therefore, 
the gauge coupling g^c runs with A^. Strictly speaking, renormalizability of NCYM is 
a subtle problem. In principle, using numerical simulations, the scaling is determined 
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nonperturbatively so that some renormalization condition is satisfied. For example, in 
|19j . D = 2 case is discussed and renormalization is performed so that the expectation 
value of the Wilson loop with the same area in physical unit is kept fixed. This result is 
equivalent to the one for the one-loop calculation. In principle, we can similarly perform 
renormalization for the case of D = A, however, it is hard with current numerical resources. 
Therefore, we rely on the one-loop calculation for this case [20]. It is known that for Non- 
Abelian gauge theory, the scaling of the gauge coupling turns out to be the same as 
that of the commutative case. On the other hand, the case of Abelian gauge theory is 
extremely different and it is known that the beta function is the same as that of Non- 
Abelian gauge theory on commutative space. That is, for NCYM, Abelian gauge theory 
is also asymptotically free as a result of the existence of non-planar diagrams. 

In order for NCYM to be well-defined, noncommutative spaces must be stable in this 
double scaling limit. However, in some cases g^c runs into a region where the space is not 
stable any more. We show that this is the case for most of bosonic models. Therefore, as 
suspected for a long time [21] , NCYM on fuzzy is not well-defined nonperturbatively 
(At least matrix model formulations discussed in this paper do not work). Here we also 
notice that D = 2 pure NCYM is only one exception that we have found in this paper. 
In other words, NCYM describes a wrong vacuum and hence noncommutative spacetime 
is not an emergent background in this case. This is not necessarily a negative conclusion 
- we can say that NCYM correctly describes gravitational instability. 

On the other hand, once supersymmetry is introduced we can expect that fuzzy spaces 
are stabilized because of the BPS nature and noncommutative super Yang-Mills theory 
(NCSYM) on fuzzy is realized. In order to formulate NCSYM on fuzzy M^, we add a 
cubic Myers term to the usual IKKT-like matrix models. One thing we notice here is the 
fact that this models themselves do not have supersymmetry but it recovers in the double 
scaling limit. 

Organization of this paper is as follows. In ^ we study bosonic matrix models to 
formulate bosonic NCYM. We firstly discuss the twisted Eguchi-Kawai model [12] and 
explain that we cannot formulate D = 4 pure NCYM [22] while we can formulate D = 2 
pure NCYM. Next, we discuss bosonic analogues of IKKT matrix models with a cubic 
Myers term and analyze the stability of solutions like fuzzy spheres. We show that we 
cannot formulate D = A and D = 2 NCYM with adjoint scalars using this formulation. 
We also demonstrate that pure D = 2 NCYM can be realized by adding a potential term 
to an adjoint scalar. In the end of this section, we comments on other scaling limits like 
commutative limit. In ^we study approximately supersymmetric matrix models with a 
cubic Myers term to formulate NCSYM and show that the approximate supersymmetry 
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stabilizes fuzzy spaces. 



2 Bosonic matrix models and bosonic NCYM 
on fuzzy 

In this section, we study bosonic matrix models and their double scaling limit which 
leads to bosonic NCYM on fuzzy R-^. In §2.1l we briefly review the twisted Eguchi-Kawai 
model (TEK) [12] and discuss the stabihty of the ground state [22l[23]. In §2 .21 we explain 
the formulation of NCYM using TEK [SI [24] and explain the double scaling limit. It turns 
out that NCYM on fuzzy cannot be realized using it [22]. In §2.3l we introduce bosonic 
analogue to IKKT matrix model with a cubic Myers term, which has fuzzy S"^ x S*^ as a 
classical solution. We show that this background is unstable in the double scaling limit. 
Discussion in this subsection applies also to other deformations with a cubic Myers term. 
In §2.41 we study other possible limits including commutative limit. 



2.1 Twisted Eguchi-Kawai model 

Twisted Eguchi-Kawai model (TEK) [12] is a unitary matrix model defined by the 
action 



Stek = -PNj2z,,uTr {U,U,UlUl) 



{21] 



where are N x N unitary matrices with the Greek indices run from 1 to D and f] is 
the inverse of the 't Hooft coupling. We mainly concentrate on the case oi D = 4. We 
comment on the case of i5 = 2 in the end of the next subsection where we discuss matrix 
formulation of NCYM on fuzzy M^. 

The phase factors Z^^^ are defined by 



n 



exp , 
In this paper, we use the skew diagonal twist which is written as 
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where L = y N corresponds to the lattice size [12] . There are other ways of twisting, but 
discussion is completely parallel and conclusion is the same as far as the double scaling 
limit which leads to NCYM is concerned. 

In the weak coupling limit {(3 oo), the path- integral is dominated by configurations 
with the minimum value of the action. This configuration 
and satisfies the 't Hooft algebra 



r(0) 



is called "twist eater" 



r r 



(2.3) 



For the skew-diagonal twist, we can easily construct a twist eater configuration by intro- 
ducing L X L "shift" matrix Sl and "clock" matrix Cl 



/ 1 
1 



Sr 



•■ 

V 1 • ■ 

These matrices satisfy 



0\ 
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0/ 



/ 1 
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g27rj-2/L 



327ri(L-l)/L 



(2.4) 



(2.5) 



and then we can construct a twist eater configuration for the above skew diagonal twist 
as 



T3 = 1l®Cl, 



T2 = Sl® 1l, 
14 = 1l ® Sl. 



(2.6) 



This twist eater configuration is nothing but fuzzy T'^ in the context of NCYM. We will 
explain the relation between fuzzy and fuzzy when we use TEK as a potential 
nonperturbative formulation of NCYM on fuzzy in the next subsection. 

In [22], it was shown, by Monte-Carlo study of TEK, that the configuration deviates 
from the and the fuzzy torus collapses in a certain range of the inverse 't Hooft coupling 
p. The upper boundary of this region scales as 



/3c ^ 0.0034iV + 0.25. 



(2.7) 



We can estimate this behavior easily and somehow roughly as follows. For simplicity, we 
assume that the fuzzy torus = F^ collapses to the identity configuration = 1^. 



The difference of energy between these configurations is 

= Stek{U^ = In) — Stek{U^j_ = F^ 



{21 
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Far from the weak coupling limit, the system has quantum fluctuations. Especially quan- 
tum fluctuations about twist-eater is known to be 0{N'^) [25]. Roughly expecting that 
the fuzzy torus collapses if the fluctuation around twist-eater conflguration exceeds the 
energy difference AS, we can estimate the critical point (3^ on which the torus begins to 
collapse as0 

(3c ~ AT, (2.9) 
which is consistent with the numerical results ( 12. 7p . 



2.2 TEK and NCYM on fuzzy 



TEK is a potential nonperturbative formulation of pure NCYM on fuzzy MP. [5| [2^. 
In order to realize the formulation, we notice that fuzzy IR-^ is realized as a tangent space 
of fuzzy . We can determine whether we can formulate the NCYM or not by analyzing 
the stability of the torus in the double scaling limit. Here we review the formulation 
of NCYM on using TEK and especially discuss the double scaling limit [22] and the 
stability of the fuzzy T^. We also comment on the case oi D = 2. 

By taking f/^ = e"^'', where a corresponds to the lattice spacing, and expanding the 
action of TEK fl2.ip . we have its continuum version as 

Stek = Yl ([^A*' - '^'-)' ' (2.10) 

up to higher order terms in a, where 

27rr) 1 

Then, by expanding the action around a classical solution (12.101) 

4°^=P/., [p„Pu]=^ef.., (2.12) 

we obtain the U{1) NCYM on fuzzy as follows. Let us deflne the "noncommutative 
coordinate" x'^ = {0~^Y^'' py. Then we have 

[x^',x''] = -i{e^^Y\ (2.13) 

This commutation relation is the same as that of coordinates on fuzzy with noncom- 
mutativity parameter 6*, and hence functions of x can be mapped to functions on fuzzy 



In this paper, we often estimate the power of N only and we use "A ^ i?" (rcsp. "A < B") to 
represent that the order of A is equal to (resp. equal to or less than) that of B. 
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More precisely, we have the following mapping rule: 

f{x)g{x) ^ f{x)'kg{x), 



(2.14) 



Tr - ^/rf 

where -k represents the noncommutative product, 

f{x)^g{x) = /(x)exp (^-'-d ,{e~'rd,^ g{x), (2.15) 

and we obtain f/(l) NCYM with coupling constant 

gj^^ = Air'^g^/Vdet 9. (2.16) 

In order to keep the noncommutative scale 6 finite, we should take the double scaling 
limit with 

a-i ~ A ~ N^/\ (2.17) 

As we have explained the identification to formulate pure NCYM using TEK, we next 
determine the double scaling limit explicitly and discuss the stability of the fuzzy T^. The 
one-loop beta function for D = 4 t/(l) NCYM is given by [20] 



Pi-ioop{gNc) = -^4^y^^c + 0{g%c)- (2.18) 
Therefore, the inverse 't Hooft coupling /5 scales as 

P ^ log A ~ log A^. (2.19) 

9nc 

Since we know that the torus collapses below the critical point which scales as (12.71) . 
we can see that the fuzzy collapses in the double scaling limit. Therefore we finally 
see that we cannot formulate D = 4 pure NCYM using TEK. 

Before closing this subsection we comments on the results for D = 2. In this case, 
(12.21) and (12. 6p are replaced by 

M =(^^^ J j . (2.20) 

and 

Ti = Cn, T2 = Sn. (2.21) 
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This corresponds to a fuzzy T^. In order to take the double scahng hmit with the noncom- 
mutative parameter 9 ~ {a^N)~^ fixed, we must scale the lattice spacing as a ~ N'^/"^. 
For D = 2, the double scaling limit is determined by numerical simulations as [19] 

/3 ~ iV. (2.22) 

In this scaling, ~ g]^^, does not run, which is consistent with one-loop beta function 
for D = 2. Since we know that fuzzy does not collapse in D = 2 TEK (This case 
is exceptional because there are no physical degrees of freedom.), we see that we can 
formulate D = 2 pure NCYM on fuzzy M^. For detailed simulations and renormalizability 
see [IS]. 



2.3 Matrix model with a cubic Myers term 

In this section, we use an bosonic analog to IKKT-type matrix model with a cubic 
Myers term. More concretely, we consider a matrix model which has S"^ x S*^ as a classical 
solution by choosing the cubic term coupling appropriately. Then we discuss stability of 
fuzzy S'^ X S'^. Although we use the specific solution, this argument itself can be applied 
to the case of other compact noncommutative manifolds like fuzzy S*^ and fuzzy CP^. 

Let us start with the d = 6 bosonic analog to IKKT model with a cubic Myers term. 
The action is written as 




(2.23) 



where is N x N hermitian matrix and the Greek indices run form 1 to 6. In the 
cubic term, f'^'^^ is the structure constant of SU{2) x SU{2) and a is a constant which 
characterizes the radii of fuzzy spheres. We choose the totally anti-symmetric tensor f^^P 
such that the only nonzero components are /^^^ = /^^^ = 1 and their permutations. 
The equation of motion for this model is 

[[A^, A,l A,] + 22«/''^M,A, = 0. (2.24) 

A classical solution called fuzzy S"^ x S"^ is given by 

4°) = aJ„ (2.25) 

^ The argument below can be parallelly applied to the case of fuzzy S'^ and fuzzy CP^ because they 
are classical solutions of (|2.23[) with f^^'^P appropriately chosen. Fuzzy S'^, however, is a classical solution 
to a bosonic matrix model with a quintic Myers term. In this case, perturbative calculation is not valid. 
In [26], it is numerically shown that >S^ is unstable in a bosonic matrix model. 
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where is a generator of SU (2) x SU (2) which satisfies 

[J,. Ju] = tr'Jp. (2.26) 

can be expressed aj^ 

^1,2,3 = ^1,2,3 ® I2S+I, ^4,5,6 = l2s+l ® ^1,2,3 (2.27) 

where J*^*) is the spin-s generator and = (2s + 1)^. In this case two fuzzy spheres have 
the same radius and square of the radius R of the fuzzy sphere is given by 

3 2 



^' = EK0 =«M^ + 1)- (2-28) 



i=l 



Expanding the matrix model fl2:23|) about (lOTj) . we obtain NCYM on fuzzy 5"^ x S"^ 
coupled to two adjoint scalars. By zooming-up the north pole, i.e. considering only states 
with J3 ~ Je ~ s, we formally obtain NCYM on fuzzy with two adjoint scalars, which 
originate from transverse directions of the fuzzy 5*^ x S*^. Because 

[Af), ] = [Af , 4°^] = 2«V3 ~ za's, (2.29) 

the noncommutativity parameter 6 is 

e r^a^sr^ a^TiV. (2.30) 

In order to keep 9 fixed, we must scale a as 

a ~ A^-l/^ (2.31) 

and therefore the momentum cutoff scales as 

A ~ as ~ iV^/^ (2.32) 

As a result, in order to take the continuum limit with 6 fixed, we have to scale as [20] 



47r2 

= -0^9 nI ~ log a ~ log AT. (2.33) 

We can easily see that fuzzy S"^ x collapses when p < A^, because energy difference 
between fuzzy S'^ x S'^ and = is of order ^^r- ~ while quantum fluctuations are 
of order A^^. Therefore, fuzzy x 5*^ collapses when we take the double scaling limit 
fl2.3ip . (12.331) and then we cannot take the continuum limit. 



We can also combine generators with different spins, but the argument does not change quahtatively. 
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This bound was derived more rigorously using Monte-Carlo simulation. Interestingly, 
this bound can also be derived through perturbative calculations of the matrix model [27] . 
Firstly, notice that eigenvalues are concentrated around the origin after the collapse of 
fuzzy sphere. This can be confirmed by numerical simulations. Then it is reasonable to 
assume that, in the perturbative analysis, such instability can be detected by considering 
only "rescaled fuzzy sphere" ^^escaicd _ ^(o) ^ const and by calculating its free energy as 
a function of the radius. At large enough coupling, there is a minimum to the free energy, 
which indicates that the background is stable. However, below some critical point, the 
minimum disappears and we can expect that the background is not stable any more. This 
is actually the case and the critical value obtained in this way agrees with numerical result 
very accurately. For details, see [27]. In the next section. We assume the validity of the 
perturbative calculation and use it to justify the matrix formulation of supersymmetric 
noncommutative Yang-Mills theory. 

Here we comment on results for the formulation of D = 2 NCYM with an adjoint 
scalar. Let us take f'^'^^ to be the structure constant e^jj^p of SU{2) where the Greek 
indices run form 1 to 3. As a classical solution of this matrix model we can obtain fuzzy 
S"^. By zooming up the north pole as we did above, we obtain the NCYM on fuzzy M^. 
In order for the noncommutativity parameter to be fixed, we have to scale the coupling 
constant for the cubic Myers term as 

a ~ iV"^/l (2.34) 

Because potential difference between the fuzzy S'^ and = is of order ~ ^ while 
one- loop fiuctuation is of order A^^, fuzzy S*^ collapses when ^ N. On the other hand, 
we can see the gauge coupling constant does not run similarly to the case of §2.2[ 
Therefore, we cannot take the continuum limit as D = 2 NCYM with an adjoint scalar. 



2.3.1 Adding potential terms for adjoint scalars 

In [28], another matrix model formulation of NCYM is introduced. This matrix model 
has fuzzy S"^ as a classical solution. In the original paper above, the commutative limit 
^ ^ oo was studied. In this section, we rather discuss the double scaling limit with 6 
fixed and see whether we can use this matrix model to formulate NCYM. 

For this model the action is given by 

S = -^,Tr I (aA, + ie,,,A,A,f + -^{N^ - 1)) | • (2.35) 
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By expanding the action about a classical solution 



A = (2.36) 

where Jj are SU{2) generators with spin s = ^^y^, NCYM on fuzzy S*^ is realized. (The 
second term in fl2.35p gives potential for adjoint scalar.) To take a continuum limit 
with fixed noncommutativity parameter, we should take large- iV limit with fixed and 
a ~ — 

However, we can easily see that this background can collapse to a point e.g. 



A^=A2 = 0, As = . (2.37) 

We can easily see the difference of tree-level potential at (12.361) and a fl2.37p is of order 
while quantum fluctuations are of order A^^ in the double scaling limit. Therefore 
we can see that the critical coupling is ^ ~ and fuzzy collapses in the limit with 

9c 

± < ± 

In [29] slightly generalized version of (12.351) . 

S = NTr {-^[X.,X,f + |^6--^X.X,X, - m^p^Xf + ^^{X!?^ , (2.38) 

was studied both numerically and perturbative and the critical point is found to be 

/ 3 \ 1/4 

By redefining the field and by identifying parameters as 

A, = g'/^N'/^X,, ^ N\ a ~ g"^N^'^p, (2.40) 

we obtain (12.351) from (12.381) up to 0(1) factors. With this identification, the scaling of 
the critical coupling becomes 



9l 



~ {N^'^p,a-^Y ~ N, (2.41) 



which agrees with the rough estimation just below (12.371) . Therefore it finally follows that 
we cannot formulate D = 2 NCYM with an adjoint scalar using (I2.35p . 

However, the generalized model ( ]2.38p has another NCYM limit. To prevent fuzzy 
sphere from collapsing in the continuum limit {g'^ fixed and a ~ we have to scale 

^ < 0(1). To realize this scaling with the redefinition of the field and identification of a 
shown in (I2.40p . we have to scale m as 

> iV^ (2.42) 
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instead of iV^. Since last two terms in (12.381) are rewritten as 

2m^ ( iV^ - 1 \ ^ 

they suppress the fluctuation perpendicular to fuzzy sphere. Therefore, an adjoint scalar, 
which corresponds to this direction decouples and we obtain D = 2 pure NCYM with the 
scahng ([232]). 

Before closing this subsection let us remark on the subtlety in the above argument. 
The bound (12.391) is obtained by calculating the free energy of the rescaled fuzzy sphere. 
However, if the value of m is extremely large, the instability (if exists) cannot be captured 
in this way, because collapse without changing the value of is more economical. If 
adjoint scalars decouple, the situation is analogous to the case of TEK. \n D = 2 TEK, 
fuzzy does not coUape. Therefore in the case of fuzzy S"^, we do not expect this kind 
of instability. On the other hand, in the case of fuzzy S"^ x S"^ or fuzzy CP^ with adjoint 
scalar potentials [31], [30], we expect this instability similarly to fuzzy in D = 4 TEK 
and hence NCYM on 7u.zy cannot be obta.ne.fl. It is desirable to check it directly 
with Monte-Carlo simulation. 



2.4 Other limits 

So far, we considered the continuum limit with the noncommutativity parameter fixed 
and showed that most of bosonic models have instability. In this subsection, we discuss 
other possible limits. For concreteness we consider D = A TEK model. 

First, let us consider the case in which the fuzzy torus does not collapse. The non- 
commutativity parameter is expressed as 

1 A2 

e ~ —= =. (2.44) 

x/iVa2 a/at 



To prevent the fuzzy torus from collapsing, the momentum cutoff must be large enough 
so that 

logA- ^ ~ A-VA^> A-^A^l /2.45) 
9 



On the other hand, to keep the volume of noncommutative space ayN ~ y N / K nonzero, 
A cannot be so large: 



A < ^N. (2.46) 



^There are models in which adjoint scalars arc dropped by hand. We expect the situation is the same 

[siiisa- 
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The only solution to the above constraints fl2.45p and fl2.46p is 



(2.47) 



up to log A corrections. In this limit, noncommutativity length goes to zero and 
spacetime volume is fixed. This limit has been studied in many references. This limit is 
of interest as an alternative to lattice gauge theory, because it might provide simpler way 
to introduce chiral fermions [53] . 

Next let us consider the case that fuzzy torus does collapse. From D-brane point of 
view, it just means D-brane collapses to lower dimensional configuration. From NCYM 
perspective this limit seems not to have a sensible continuum limit because there is no 
extended direction. In [3l] slightly different model with two commutative and two non- 
commutative dimensions has been studied numerically. In that case two noncommutative 
dimensions collapse similarly to our case, but numerical results suggest that there is a 
continuum limit with two commutative noncompact directions and two compact, finite 
size "noncommutative" directions. Such models would be interesting as a toy model for 
compactification mechanism in matrix models. 

3 Supersymmetric matrix model and noncommuta- 
tive super Yang-Mills 

In the previous section, we have discussed various matrix model formulations of bosonic 
NCYM. In this section we explain the formulation of noncommutative super Yang-Mills 
(NCSYM). For this purpose we introduce matrix models with an approximate supersym- 
metry and perturbatively discuss stability of noncommutative spaces. 

Let us consider the IKKT-like matrix model [3] with a cubic Myers term 



where and ip are bosonic and fermionic Hermitian SU {N) matrices, Greek indices run 
from 1 io d {d = 4,6,10), ip has a spinor index and F'^ is the SO{d) Gamma matrix. 
j/iup ^YiQ structure constant of a Lie group whose rank r is less than d. Except for the 
cubic Myers term, we can obtain this action from D-dimensional A/" = 1 SU{N) super 
Yang-Mills by dimensional reduction. 

Since numerical simulations for these matrix models are difficult, except for d = 4 case 
[35] due to the notorious sign problem, it is difficult to discuss stability of backgrounds 





(3.48) 
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nonperturbatively. Hence, we provide only perturbative arguments, which works perfectly 
well for bosonic models. The perturbative argument is carried out similarly to the case 
of bosonic analog of IKKT-like matrix model with a cubic Myers term. 

Although this model is not supersymmetricjf], the noncommutative background can 
be stabilized for any value of a. At large a it is stable because potential barrier is very 
high and, furthermore, fluctuations are suppressed due to approximate supersymmetry. 
At small a, it can be stabilized since this model is almost supersymmetric (at a = the 
supersymmetry recovers) and this background is almost BPS. 

As a concrete example, let us take d = 10 and ] to be the structure constant of 
SU{2) X SU{2). Fuzzy S'^ x S'^ ^7n\ is one of the classical solutions for it. (Indeed, there 
is a subtlety for this background. We will discuss it in §3.1[ ) Quantum corrections to 
this background is calculated in [3H]. Here we consider the deformation in radial direction 
only as we have explained in §2.3[ Up to the leading order of 1/A^, the tree level action 
Vtree and the one-loop correction V\^ioop for rescaled fuzzy sphere = (1 + e)aJf^ are 
calculated as 

ri_,„^ = ]V-21og2. (2 + ^^^^) . (3.49) 

(See Appendix |A] for derivation.) If we scale a ~ N~'^^^, we have 

Ttree ~ -Nq'^ T^_loop ~ A^. (3.50) 

The matrix model we are considering here is expected to realize D = 4 J\f = 4 NCSYM in 
the continuum limit and then the coupling g does not run at one-loop level in this limit. 

From (13.501) . the one-loop correction is smaller than tree- level action provided that 
is sufficiently small. We also notice that n-loop effect is 0{N{g'^/{a^N))^~^) = Ng'^^"'~^^ 
as a result of approximate SUSY and higher loop effects are negligible in this case [38] . 
We therefore see that the classical minimum e = survives after taking into account 
quantum corrections and we can expect that the fuzzy S"^ x S"^ does not collapse. On 
the other hand, at strong coupling it might collapse. To overcome this difficulty, it is 
probably useful to consider a supersymmetric deformation in [3^ HO] . 

In the above construction using the matrix model with the cubic Myers term, only 
extended supersymmetry can be realized. In order to construct J\f = 1 NCSYM, super- 
symmetric generalization of TEK would be necessary. 

^ d = r = 3 model with the cubic term is supersymmetric j36j . However, in this case the finitencss of 
partition function is not known for generic [37] . 
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3.1 Subtlety for x 5^ case 

In this subsection we discuss a subtlety for fuzzy S*^ x background. 

Because fuzzy S"^ has smaller free energy, fuzzy S*^ x S"^ is not stable; one of the 5*^ 
can shrink, while the other expands [H]. We notice that SU{2) x SU{2) is preserved in 
this proces and that this instability cannot be read off from the one-loop effective action 
f l3.49p . To avoid this instability, we should use four-dimensional fuzzy manifolds with 
higher symmetry, e.g. the fuzzy CP^. CP^ can be stable since the symmetry must be 
broken during the transition to 5^. The effective action does not change qualitatively [12] 
and we can realize D = 4 Af = 4 NCSYM using CP^. 

It is difficult to realize D = A M = 2 NCSYM using matrix model formulation because 
of the instability of fuzzy S"^ x S"^. Naively, if we add a cubic Myers term to d = Q 
supersymmetric matrix model as above, four-dimensional M = 2 NCSYM is expected 
to be realized in the continuum limit. In this case, the coupling runs as g""^ ~ logA^, 
and hence the background is stable. However, to realize four- dimensional compact fuzzy 
space with 6 matrices in this model, we need to use S*^ x S"^. It is necessary to fix the 
radii somehow, for example by quenching the background or adding a small potential 
term to the adjoint scalars, while keeping the continuum theory unchanged. Instead of 
fuzzy S'^ X S"^ the fuzzy S'^ might be useful. To make fuzzy S*^ a classical solution, we 
have to add quintic Myers term. However, it is difficult to discuss the stability because 
perturbative calculation is not valid. 

4 Conclusions and discussions 

In this paper, we studied the stability of noncommutative spaces in several matrix 
models and discussed whether or not they provide nonperturbative formulation of non- 
commutative Yang-Mills theory (NCYM). It turns out that most of matrix model formu- 
lations of bosonic NCYM on fuzzy do not work. The only exception we found is -D = 2 
pure NCYM. In the context of D-branes dynamics, those not realized correspond to false 
vacua. This might be a negative conclusion if one regards NCYM itself as a UV complete 
theory. However, as an effective description for a D-brane system, these bosonic NCYM 
correctly reproduce the instability of the system. According to [21], large one- loop cor- 
rection to free energy, which leads instability of NCYM, is due to UV/IR mixing. Hence 
by eliminating UV/IR mixing somehow, we expect that NCYM be stabilized. 

On the other hand, noncommutative super Yang-Mills(NCSYM) on fuzzy with ex- 
tended supersymmetry can be formulated using a supersymmetric matrix model deformed 
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by a cubic Myers term. At least, as we have see above, D = A JV = 4 NCSYM in weak 
coupling is realized using this formulation. Also in certain non-supersymmetric model 
with adjoint fermions, Z^v symmetry is not broken [43]. Then combining it with twist 
prescription a certain non-supersymmetric NCYM will be obtained. 

Here we comment on the formulations of NCSYM at finite temperature. For this pur- 
pose, we consider supersymmetric matrix quantum mechanicj^ with Euclidean time direc- 
tion compactified and antiperiodic boundary condition for fermionic variables imposed. At 
high temperature, fermionic modes decouples and the theory becomes essentially bosonic. 
Therefore, we can expect that noncompact fuzzy space cannot be constructed in the high 
temperature limit. Whether NCSYM at nonzero temperature exists or not is a subtle 
problem and numerical simulation along the line of [HI HS] will be necessary. 

Though we have discussed matrix models formulation only in this paper, there is 
another candidate for nonperturbative formulation of bosonic NCYM [IS]. However, 
it seems to share the same problem with the matrix model formulation considered in 
this paper. In [18] NCYM is mapped to a lattice gauge theory with twisted boundary 
condition. In the continuum limit with noncommutativity parameter fixed, however, 
corresponding lattice gauge theory goes to zero volume and essentially reduces to the 
TEK model (see Appendix iBl) 

Of course, the pathology discussed above does not prevent us from nonperturbative 
formulations of non-gauge theories on noncommutative spaces using matrix models. For 
example, scalar field theories are well defined and we can numerically analyze them using 
matrix model formulations [19][50]. We also notice, as explained in §2.41 we can take the 
"commutative" limit of NCYM, in which the noncommutativity length goes to zero. 
Therefore, one may still regard NCYM as an alternative to the lattice construction for 
gauge theories on commutative spaces. 

In the end, we comment on the recent progress in TEK and its relation to matrix 
formulation of bosonic NCYM. Since collapse of fuzzy sphere in TEK model is nothing 
but the breakdown of Z^r symmetry (original motivation for TEK is to keep this symmetry 
unbroken), construction for bosonic NCYM is tightly related to a modification of Eguchi- 
Kawai model such that Z^r does not break and large- reduction works. Historically 
two options have been studied. One is TEK, which works fine ai D = 2 but turns 
out to fail at Z) = 4. Another one is the quenched Eguchi-Kawai model (OEK)[TT]. 
in which commutative and extended background is "quenched" by hand. Naively by 

^ Monte Carlo simulation for supersymmetric matrix quantum mechanics without cubic term has been 
performed recently |441 145j. and incorporation of a cubic term |46| will be straightforward. Thermody- 
namical property of fuzzy sphere in bosonic model is studied in [47]. 
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combining twist and quench prescriptions, i.e. by fixing noncommutative background by 
hand, NCYM seems to be reahzed. However, it does not seem to work. Indeed, recently 
it was argued that QEK does not work due to the following reason [51]. In QEK, unitary 
link variables t/^'s are constrained to be V^e^^f^V^, where = diag(p^, ■ ■ ■ ,p^) is fixed, 
V^'s are unitary matrices and p^'s are distributed uniformly in M^. Naively one expects 
V^'s fluctuate around Iat and, therefore, Z^r is not broken. However, what actually 
happens is that V^'s become certain permutation matrices, so that quenched momenta 
are "locked" [51] and free energy becomes smaller. Intuitively, this result implies, even 
if the background is quenched by hand, can get a nontrivial VEV and an essentially 
different background emerges. 

The same can take place also when we quench the noncommutative background. Such 
a subtlety does not exist in supersymmetric case, and D = A M = 2 NCSYM would be 
realized by quenching fizzy S"^ x 5*^ background. 

Recently a new deformation to Eguchi-Kawai model was proposed in [52] . They added 
potential terms for Wilson lines to prevent Z^r from breakdown and argued that the 
additional terms do not contribute in the continuum limit. If it really works, by combining 
this method with the twist prescription, bosonic NCYM might be realized. Then, it would 
be interesting to understand the meaning of the deformation in the context of D-brane 
dynamics. 

In this paper, we assumed the running of the coupling constant is determined by one- 
loop beta function when we discuss the case of D = 4. However, renormalizability of the 
NCYM is of course controversial. It will be better if we can determine the running more 
rigorously, for example, by calculating correlation functions using numerical simulations. 
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A Derivation of one-loop effective action in super- 
symmetric matrix model 

Let us expand the action 

S = \Tr (-\[A,, A,]' + '^ar^A.A^A, - ^^r^[A^, ifj]) (A.51) 
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about the rescaled fuzzy sphere 



= (1 + e)aJ^. (A.52) 

At tree level, we have 



^tree = ^Nsis + 1) |(1 + e)' - "(l + j ^ r__ |(i + - -(1 + } , (A.53) 
where = (2s + 1)^. Then, the one-loop effective action is [3] 



2 ° V 1 + e (adjy J 4 [ V 2 (adJ) 

(A.54) 

To leading order in A^, we have 



\ 1 
2 + Tr- 



1 + e) V (adjy 



= A^-21og2-(^2 + ^^^), (A.55) 
where we have used [38] 

= Tr ^ J2J + W+1L . 2Ariog2 

(A.56) 



B Lattice formulation 

Lattice regularization [IS] relates commutative U (N) lattice gauge theory on twisted 
torus to a "lattice regularization" of U{1) NCYM on periodic fuzzy torus. Basically this 
relation IS cLS cL result of the fact that the Morita equivalence holds at lattice level. 

For simplicity, we consider the D = 4 U{N) gauge theory on a rectangular four-torus 
with period L. The action is 

S = [Uf.ix)U,ix + ail)U^{x + auyU^ixY] , (B.57) 



18 



where [/^ are unitary matrices which correspond to U{N) gauge fields. They satisfy 
twisted boundary condition 



[B.58) 



where Ti, are twist-eaters appeared in §2.11 

We now introduce a map A(x) between lattice fields U^{x) and operators IJ^ as 



(B.59) 



where the mapping function A{x) is defined as 



A(x) 



m'-eZ/n \i=l 



(B.60) 



where ka is a momentum ka = 2'Kma/l and n is a integer n = l/a. 

In order to relate operators to noncommutative U{1) gauge fields, we now introduce 
another mapping function A'(x') defined as 



A'ix' 



N 



(BM] 



where /' = IVN, k'^ = 27rma//', n' = I' /a and 



e 



ab 



/ \ 

-0 



V 0-00/ 







[B.62) 



ab 



We have used primed quantities to represent those on a lattice corresponding to A'. 
This A'(x') maps the noncommutative lattice fields to operators whose dimensionless 
noncommutativity parameters is 0. Because of the twist boundary condition of U^{x) 
the operator have another expansion using A'(x'), 



(B.63) 



where U'^{x') are noncommutative f/(l) gauge fields which live in periodic torus whose 
size is /' and the dimensionless noncommutativity parameter is 0. 
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Now we gain a map from U{N) gauge fields U^lx) on a twisted commutative torus to 
the noncommutative U{1) gauge fields U'^{x') on a periodic fuzzy torus. Indeed the action 
( ]B.57p is rewritten in terms of U'^{x') as 



^ = -4 E E * UUx' + afi) ^ U'^ix' + auY ^ UUx'Y] , (B.64) 



Q 



where 



/ = ^i- (B.65) 



Dimensionful noncommutativity parameter, which appears in commutators of coordi- 
nates is given by 

= Q.^ = L^, (B.66) 

Now let us consider the limit which leads to fuzzy with finite value of 6. To fix 6, 
we have to take 

/ ~ iV-l/^ (B.67) 

that is, we have to take infinitely small twisted torus and the model essentially reduces to 
TEK. Therefore it is plausible that the center symmetry U{1) breaks down. This means 
that the fuzzy torus collapses and we cannot realize fuzzy M"^ which is expected to appear 
as a tangent space of the torus. 
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